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The first measurements of ligfit antinucleus production in Au+Au collisions at RHIC are reported. 
The observed production rates for d and '^He are much larger than in lower energy nucleus- nucleus 
collisions. A coalescence model analysis of the yields indicates that there is little or no increase 
in the antinucleon freeze-out volume compared to collisions at SPS energy. These analyses also 
indicate that the "^He freeze-out volume is smaller than the d freeze-out volume. 



The Relativistic Heavy-Ion Collider (RHIC) at 
Brookhavcii National Laboratory has recently begun op- 
eration with Au beams at ^/spfj\j — 130 GeV and extends 
the available center of mass energy in nucleus-nucleus col- 
lisions by nearly a factor of 8 over CERN SPS collisions 
at y/s 



P/A through a coalescence parameter Ba, 



NN 



17 GeV. First measurements from RHIC 
indicate an increase of at least 70% in the charged mul- 
tiplicity for central collisions compared to previous mea- 
surements [|l| . Measurements of the antiproton to proton 
ratio at mid-rapidity indicate that the central collision 
region is approaching the net-baryon free limit. Such a 
system with large multiplicity and small net-baryon den- 
sity is well suited for the production of light antinuclei. 
In this letter, we report the first measurements of d and 
"^He production at RHIC. 

At RHIC energies, production of antinuclei is possi- 
ble via two mechanisms. The first mechanism is di- 
rect production of nucleus- antinucleus pairs in elemen- 
tary nucleon-nucleon or parton-parton interactions. The 
RHIC center-of-mass energy is well above the threshold 
for such processes. Due to their small binding energies, 
nuclei or antinuclei produced via early direct production 
are likely to be dissociated in the medium before escap- 



The second, and presumably dominant, mechanism for 
antinucleus production is via final-state coalescence ^, 
p. In this picture, produced antinucleons merge to form 
light antinuclear clusters during the final stages of kinetic 
freeze-out. The measured yield of nuclei or antinuclei 
with nucleon number A and momentum P is related to 
the primordial nucleon invariant yield at momentum p = 



E 



(PP 



Ba{E- 



(1) 



Equation ^ requires that antineutrons and antiprotons 
are produced with identical momentum spectra. 

Previous studies of smaller collision systems have noted 
that the measured coalescence parameter Ba can be di- 
rectly predicted from the nuclear wave function of the 
produced (anti)nucleus When going to higher ener- 
gies or larger collision systems, however, the measured co- 
alescence parameter is lower than that measured in small 
systems. This can be understood by noting that once the 
collision region is larger than the intrinsic size of the pro- 
duced (anti)nucleus, (anti)nucleons of equal velocity are 
not always in close proximity and hence do not always 
form a bound state j6j. In this sense, the coalescence pa- 
rameter can be used to infer the space-time geometry of 
the system [Q. Measurements of light nuclei and antin- 
uclei are thus analogous to two-particle Hanbury-Brown 
Twiss correlations (HBT) in that they measure "homo- 
geneity lengths" of the system at kinetic freeze-out Q . 

The measurements were made using the STAR detec- 
tor 1^ . The main tracking detector is a cylindrical Time 
Projection Chamber (TPC), which resides in a solenoidal 
magnet that was operated with a field strength of 0.25 
T for the data reported here. The TPC tracks and 
identifies most charged particles produced in the cen- 
tral pseudo-rapidity region (—1.8 < rj < 1.8) with nearly 
full azimuthal coverage. Events are selected on the basis 
of coincidence of spectator neutron signals in two Zero- 
Degree calorimeters located ±18.25 m from the nominal 
interaction region. Central events are selected using a 
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Central Trigger Barrel (CTB) that measures the charged 
particle multiplicity with full azimuthal coverage in the 
pseudo-rapidity region — 1 < r; < 1. This analysis fo- 
cuses on semi-central events, where the centrality cor- 
responds to roughly the most central 18% of the mea- 
sured minimum-bias multiplicity distribution. The anal- 
ysis uses « 600,000 events where the interaction vertex 
is within the range covered by the TPC (-200 < z < 200 
cm). 

Particle identification is done by measuring the average 
ionization energy loss {dE/dx) for each track. Studies of 
the STAR electronics response show no evidence for sat- 
uration below 30 times minimum ionizing. For tracks 
with sufficient transverse momentum to leave the TPC, 
the path length exceeds 1.4 m. For the tracks used in 
this analysis, the dE/dx resolution is ~11%. For each 
track, up to 45 ionization space-point samples are taken 
along the path through the TPC. Space-points are found 
by identifying local maxima of the ADC distribution. 
Merged ionization clusters, where multiple tracks con- 
tribute, are identified by looking for multi-peaked struc- 
ture in the ADC distribution. For the current analy- 
sis of relatively rare particles, it is necessary to impose 
tight cuts to eliminate background tracks with improp- 
erly measured dE/dx. We require a track to have at least 
35 of the 45 possible space-points. For central events, 
cluster merging is quite common and can lead to prob- 
lems with the particle identification. To avoid these prob- 
lems, we eliminate potentially merged clusters from the 
sample used to calculate the dE/dx. For the final sam- 
ple, we require that no more than 30% of the measured 
space-points come from potentially merged clusters. To 
avoid the Landau tails in the dE/dx spectrum, we use a 
truncated mean of the lowest 70% of the measured dE/dx 
samples. Figure |l| shows the measured truncated mean 
dE/dx versus the magnetic rigidity for the negatively- 
charged tracks considered in this analysis. 

Figure [l] also shows the Bethe-Bloch expectation for 
d , t and ^He. There is a clear d band below rigid- 
ity « 1 GeV/c. This analysis uses only the kinematic 
region of good d particle identification and efficiency 
(0.5 <PT < 0.8 GeV/c and rapidity |yj_< 0.3). We ob- 
serve 14 counts clustered around the ^He expectation in 
the kinematic range 1.0 < pt < 5.0 GeV/c and \y\ < 0.8. 
Note that we plot the rigidity, so the momentum of the 
"^He candidates is twice as large. No clear t band is 
observed, but if one assumes that i and ^He are pro- 
duced in similar numbers and with similar momentum 
distributions we would expect the bulk of the t to have a 
higher rigidity where our dE/dx resolution is inadequate 
for their identification. 

To extract the d yield, we construct a quantity Z = 
log{[dE /dx]/ I^{p)), where /j(p) is the expected ioniza- 
tion for a J of momentum p. For a pure sample of J, this 
quantity should be well described by a Gaussian centered 
at zero. In the insert of Fig. |^, we plot the Z distribution 
for one transverse momentum bin. We see a Gaussian d 
signal superimposed upon a background due to the tail 



of the p distribution. We parameterize the p background 
in the tail region as an exponential, and fit the resulting 
distribution to a Gaussian d signal -I- exponential p tail 
hypothesis. In the insert of Fig. |^, we also show (by the 
curve) our exponential-f Gaussian fit. In the d kinematic 
region considered, the signal to background ratio ranges 
from 30 in the lowest pr bin to 3 in the highest pT bin. 
We have performed a similar analysis of the ^He Z dis- 
tribution, and estimate the total background to be less 
than 0.5 counts. For extracting yields, we assume that 
our 14 observed '^He are background free. 

To evaluate the efRciency, we use GEANT and a TPC 
response simulator to create raw pixel level simulated 
tracks which we then embed into real events. The em- 
bedding is crucial for this analysis since it allows us to 
estimate the effects of cluster merging on our efficiency. 
No data on d and ^He interactions in material exist in 
the literature, and these antinuclei are not incorporated 
into GEANT. Instead we use d and "^He simulations in 
GEANT to understand our acceptance and tracking ef- 
ficiency. We then add a correction for the estimated an- 
nihilation in the detector, where we assume that the d 
annihilation cross-section is 1.4 times the p annihilation 
cross-section, and that the "^He annihilation cross-section 
is twice the p annihilation cross-section. The p annihi- 
lation correction was discussed in a previous publication 

, and the cross-section scaling relations are taken from 
Ref. ||l^. Final calculated efficiencies are in the range 
0.2-0.5. This is much lower than the typical STAR effi- 
ciency for charged particle tracking (0.8-0.9). The differ- 
ence is due entirely to the restrictive track cuts used in 
the current analysis to eliminate backgrounds. 

Systematic errors were estimated by varying the cuts 
used in the analysis. These variations include changing 
the number of hits for a valid track, changing the allowed 
region of vertex locations, changing the assumed annihi- 
lation cross-sections, and changing the Z range used for 
the signal-|-background fit. We estimate the maximum 
systematic error on the invariant yields to be around 15%. 
We also assume that the errors on the individual yields 
are largely correlated. This causes the systematic errors 
to partially cancel when forming coalescence ratios. 

We extract d invariant yields in three transverse mo- 
mentum bins, where each bin has w 100 entries. The 
extracted yields are listed in Table |[ Comparing these 
yields to lower energies, there is a factor of ~ 50 increase 
in the d production rate in going from a/s^^ = 17 GeV 
|ll| to \/s[^]^ — 130 GeV, and an even more dramatic 
factor of w 60,000 increase in the d production rate rel- 
ative to AGS energy {y/sj^j^ = 4.9 GeV) 

The mean transverse momentum of the observed ^He 
sample is «2.4 GeV/c. We extract an invariant yield 
per event evaluated at the mean pT of [8.4 ± 2.3{stat.) ± 
1.3(s?;s.)] X lO^^GeV^^c"^. We have assumed an expo- 
nential transverse mass distribution to calculate a cross- 
section weighted average efficiency in the STAR accep- 
tance. NA52 has reported two '^He in minimum-bias 
Pb-HPb collisions at the CERN SPS Our invariant 
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yield is higher, but quantitative comparison cannot be 
made because of the different centrahties. 

Although only 14 counts were observed, our large 
kinematic coverage for "^He allows us to estimate the 
dN/dy and inverse slope T. To do this, we have cal- 
culated the expected yield as a function of y and px 
using efficiency calculations from embedded data and 
assuming an exponential transverse mass distribution. 
We minimize the negative log-likelihood over the entire 
STAR acceptance taking into account phase-space cells 
with no observed counts. Wc extract '^He dN/dy — 
[5.1 ± 1.7{stat.) ± 0.8(sys.)] x 10^^ and an inverse slope 
r = 0.70±0.25(stat.) GeV. 

STAR has measured invariant yields for p in a simi- 
lar centrality range |l^. These results can be combined 
with the invariant yields presented in this paper to calcu- 
late coalescence factors using Equation |l]. In the coales- 
cence picture, only antinucleons produced directly from 
the source are available to form light antinuclei. Hence, 
the p yields in the coalescence ratio have been corrected 
for antihyperon feeddown. We use the RQMD model 
1^ and a detector simulator to evaluate the probability 
of incorrectly assigning a weak-decay produced p to the 
primary vertex, and find that about 45±5(sys.)% of our p 
sample comes from antihyperon feeddown. This fraction 
is consistent with preliminary STAR measurements of the 
A/p ratio. Table ^ lists the total p invariant yields along 
with the estimated correction for antihyperon feeddown. 

For the top 18% most central collisions, wc find {B2) = 
[4.5 ± 0.3{stat.) ± l.Oisys.)] x IQ-* GeVVc^ in the d 
kinematic region 0.5 < pr < 0.8 GeV/c and |?;| < 0.3. 
In the top panel of Fig. || we compare this result to 
previous measurements at lower energies. Here we plot 
the results for both d and d. In pA collisions, B2 is es- 
sentially independent of the collision energy. In central 
nucleus-nucleus collisions, however, the coalescence fac- 
tor B2 decreases as the collision energy increases from 
Bevalac to AGS to SPS. The STAR result shows that 
there is no similar decrease in B2 from y/s^j^; = 17 GeV 
to -x/i^jv = GeV. Comparing the STAR result to the 
average of the two d resuhs at the SPS ^ |l|| , we obtain 
B2{SPS)/B2{RHIC) = l.l±Q.l{stat.). 

For the top 18% most central collisions, we find (-83) = 
[2.1 ± Q.&{stat.) ± 0.6(sys.)] x lO"'^ GeV^/c^ in the ^TU 
kinematic region 1.0 < px < 5.0 GeV/c and \y\ < 0.8. 
Once again, we compare to collisions at lower energies in 
the bottom panel of Fig. ||. The qualitative trend for B^ 
is very similar to S2. For pA collisions, the coalescence 
factor is independent of energy. For AA collisions, the 
coalescence factor decreases with increasing collision en- 
ergy. The statistics of the ^He measurement at the SPS 
preclude a quantitative comparison. If we compare to 
the average of '^He and '^He at the SPS [Q, we obtain 
B3{SPS)/Bz{RHIC) = 2,A±l.h{stat.). 

Several prescriptions have been proposed for relating 
the coalescence parameters to a geometrical source size 
[^, 0, ^, |l|. For these models, the coalescence param- 
eter scales with the volume as Ba oc l/y*^"^"^^ in the 



limit of an (anti)nucleon volume much larger than the 
intrinsic size of the produced (anti)nucleus. Using this 
simple expression, and the measured coalescence param- 
eter ratios, we see that V^{RHIC) = (l.liO.l) V^{SPS) 
and VsjY^iRHIC) = (1.8 ±0.4) V,jy^{SPS). Both mea- 
surements indicate no large increase of the antinucleon 
freeze-out volume when going from \/S]\[j^ = 17 GeV 
to \/i^7v = 130 GeV. STAR has also measured source 
sizes using tt^tt" interferometry p^ . If we construct a 
quantity proportional to the volume, T/^^ cx R^R^, and 
compare to the published SPS data [|2^, we estimate 
V^~^-{RHIC) = (1.8 ±0.7) V^-^~{SPS). AU three 
available measurements indicate only a slight increase 
in volume compared to lower energy collisions. Caution 
should be exercised, however, when making quantitative 
comparisons between the volumes measured via coales- 
cence and the volumes measured via HBT since it is not 
clear that the freeze-out space-time geometry for pions 
and antinucleons should be the same. 

We can also make quantitive estimates of the freeze-out 
geometry within the context of a particular coalescence 
model and ask whether the d and ^He sources are the 
same. To do this, we use a simple thermal model [0. 
This model assumes that antinucleons and antinuclei are 
in chemical and thermal equilibrium within a volume V . 
From this model, we extract = 6700 ± 500(stat.) fm'^, 
and V^-jj-^ = 3800±500(stat.) fm^. This discrepancy indi- 
cates that the thermal model assumptions are not valid 
in the production of light antinuclei. The '^He freeze- 
out from a smaller volume and at a presumably earlier 
time compared to d . This trend of decreasing source size 
with increasing nucleon number has been observed before 
in the production of light nuclei [Q, |l^ . The coalescence 
picture of light antinucleus production would predict that 
the probability for producing an antinucleus with mass 
A is proportional to the A^^ power of the local antinu- 
cleon density. If the antinucleon source is not of uniform 
density, one would expect the different mass antinuclei 
to measure different source sizes, and this is indeed what 
we observe. We have applied other coalescence models 
with different assumptions. The Sato and Yazaki model 
[^ indicates a similar trend as the thermal model, with 
Vj/VS)-^ — 2.2 ± 0.3, while the Scheibl and Heinz model 
[p|, which can be calculated assuming a Gaussian antin- 
ucleon density profile and explicitly includes the effects 
of radial flow, gives V^/Vsjj-^ = 0.9 ± 0.1. In the Scheibl 
and Heinz model, an equivalent effective volume, as indi- 
cated by the data, would imply a larger total volume for 
d compared to ^He. 

In summary, we have made the first measurements of 
the production of light antinuclei {d and ^He) in Au+Au 
collisions at \fsj\jj^ = 130 GeV. A large enhancement in 
production rate is observed compared to lower energies. 
We have combined the measured yields with measure- 
ments of p production to extract coalescence parameters 
B2 and -B3. Quantitative comparisons to SPS results in- 
dicate little or no increase of the antinucleon freeze-out 
volume. We also find that the ^He are produced from a 
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smaller volume than the d . 
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Pt(GcV/c) 


E^(GeV-2c3) 


P E^(GeV-^c3) 


w.d. correc 


d 


0.55 


[2.47 ± 0.26] X 10-3 


4.20 ±0.12 


0.56 


0.65 


[1.87 ±0.19] X 10-3 


4.00 ±0.10 


0.53 


0.75 


[1.93 ± 0.20] X 10-3 


3.82 ± 0.09 


0.52 




2.4 


[8.4 ± 2.3] X 10-'^ 


2.63 ± 0.04 


0.61 



TABLE I: Measured invariant yields of aiitiimclci. The errors 
quoted are statistical only. Systematic errors arc estimated to 
be 15%. Also listed are p invariant yields at the same velocity, 
and the weak-decay correction to the p yield estimated from 
RQMD. 
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FIG. 1: Ionization energy loss (dE/dx) versus rigidity 
(I momentum/nuclear charge unitsj) for negative tracks. The 
Tf- and K~ bands have been suppressed. Also plotted are 
the Bethe-Bloch expectations for d, t and "^He. Inserted is 
a projection of the Z variable (see text) for one transverse 
momentum bin (0.6 < pr < 0.7 GeV/c). 



10 

o 

CM 
> 

0) 

O 



Bevalac SPS 
AGS 

pA Average 



m 



10 



10 



o 



o 



8 



STAR 



10 

10-' 



> 

O 10 



pA Average 



DQ 



10 



10 



10 



10 



O 

o 



10 



10^ 

^NN [GeV] 



FIG. 2: Coalescence parameters B2 and -63 excitation func- 
tions for nuclei (hollow markers) and antinuclei (solid mark- 
ers) in semi-central Au-|-Au or Pb-|-Pb collisions |l^, |l^, 
PI 1^. The errors on the STAR data points are 
statistical (narrow bars) and systematic (wide bars). 



